Control of Raman lasing in the non-impulsive regime 
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We explore coherent control of stimulated Raman scattering in the non-impulsive regime. Optical 
pulse shaping of the coherent pump field leads to control over the stimulated Raman output. A 
model of the control mechanism is investigated. 
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Since the advent of pulsed lasers, there has been ex- 
tensive experimental and theoretical work on stimulated 
Raman scattering. More recently there has been a resur- 
gence of interest in impulsive Raman scattering^, 13,0,0 
due in part to the development of the field of learn- 
ing coherent control 0, I^JSi along with advances in 
ultrafast laser technology!^ and programmable pulse 
shaping [glliol. In the impulsive regime, the laser band- 
width is large compared to the frequency of the Raman 
mode (typically a molecular vibration). In the time do- 
main the pulse duration is short compared to the vibra- 
tional period, and one can shape the optical pulse to 
drive molecular mode(s) on resonance|llj. Alternatively, 
if the laser bandwidth contains photon pairs separated 
by the Stokes frequency, then the Raman gain can be 
seeded by the pump laser. Control of the Raman gain 
can be achieved by appropriately phasing colors in the 
pump pulse. 

For the case of two-photon atomic absorption, this con- 
trol mechanism has been described as shaping of the non- 
linear power spectrum of the driving field |l2L[l3l ll4l|. The 
idea has since been extended to multi-photon absorption 
in molecules 0, , as well as vibrational Raman exci- 
tation in multi-mode molecular systemslllj. Here again 
analysis leads to an explanation of the control via the 
nonlinear power spectrum of the pump pulse. 

In the non-impulsive regime, the laser bandwidth is 
small compared to the frequency of the mode. Since there 
are no photon pairs separated by Stokes frequency in the 
pump pulse, the Stokes field must build up from spon- 
taneous Raman scattering. This produces a Stokes field 
whose phase is random and cannot be controlled. Control 
may still be possible over the stimulated output spectrum 
or the final state populations in the presence of multiple 
Raman modes. Non-impulsive control of the stimulated 
Raman spectrum in liquid methanol has been reported 
in experiments that used lear ning algorithms to discover 
the optimal driving field [13, llSl ■ Here we demonstrate 
a possible mechanism for this control. 

The details of our adaptive learning technique and 
laser system have been described previously [ig. In 
learning control experiments, the physical system is used 
to find the optimal driving field (pulse shapes) without 
prior knowledge of the Hamiltonian ^ . These solutions 
can be analyzed later to learn about the underlying quan- 
tum dynamics. Briefly, our experiments use a shaped, ul- 



trafast Ti:Sapphire laser system as the excitation source. 
The laser pulses are shaped in an acousto-optic Fourier 
filter interfaced with a computer. Our spectral band- 
width (4 — 5 THz) and pulse shaping characteristics pro- 
vide temporal control over the pulses ranging from 100 fs 
to 5 ps in duration. Our adaptive learning algorithm de- 
termines the pulse shapes by optimizing a feedback signal 
derived from the Raman spectrum. 
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FIG. 1: Resulting Raman spectra obtained after optimization 
with the feedback algorithm. Inset shows a Husimi plot of 
the optimal "double-blob" pulse shape that led to S-mode 
excitation. 

The system under investigation in these experiments 
is liquid phase methanol (CH3OH). We focus on the two 
C-H stretch vibrational modes labelled S (2834 cm~^) 
and A (2946 cm~^) 19]. The experiments are performed 
in a 10 cm liquid cell, and the forward scattered Raman 
spectrum is collected and fed back to the algorithm. The 
learning algorithm is able to channel the gain into either 
of the two Raman modes. Although the algorithm finds 
a variety of pulse shapes that produce the desired effect, 
one class of solutions stands out. Figure ^ (inset) shows 
a representative pulse shape of this class of solutions 
in the form of an optical Husimi distribution '2^. The 
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Husimi plot is a two-dimension convolution of the Wigner 
function (20l Ell of the pulse shape that approximates 
the intensity distribution in both time and frequency. 
The optimal pulse shapes found by the algorithm con- 
sist of a pair of nearly transform-limited "blobs" that are 
separated by approximately 3.3 THz. This energy split- 
ting matches the S-A mode spacing. Figure ^ shows the 
resulting Raman spectrum after excitation by the pulse 
shown in the inset of Fig. ^ The pulse produces strong 
selective excitation of the S vibrational mode. 
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FIG. 2: Calculated Wigner plots for double-blob pulse used 
in the experiment when the phase offset between the blobs is 
equal to tt. The lower plot is the Husimi distribution when 
the phase offset is zero. 

These two-blob solutions suggest a quasi-impulsive 
model for the interaction, where the pump field directs 
the Raman gain through a coherent coupling between the 
two modes This coherence modulates the intensity 
envelope of the pulse. The modulation is apparent in the 
Wigner distribution, but has only a minimal effect on the 
Husimi distribution (Fig. 2). 

In order to study the effect of the phase offset, we ex- 
cite the molecules with the original pulse shape shown in 
Fig. 13 We then collect the stimulated Raman spectra 
as the phase offset between the two blobs is varied from 
to TT. Figure O shows the stimulated Raman output 
spectrum for several different phase offsets. As the off- 
set is varied, the Raman gain oscillates between the two 
different Stokes modes. As is to be expected, this single 
parameter control is not as effective at mode selection 
as the full learning algorithm (see Fig. Nevertheless, 
control is achieved by varying only the phase offset. Care 
was taken to adjust the laser intensity to avoid saturation 
of the Raman gain, since this leads to significant changes 
in the observed control. 

A simple model that predicts the observed effect can 
be constructed by expanding on previous work on single 
mode Raman scattering. The theory governing stimu- 
lated Raman gain of a single active mode under exci- 
tation by an off-resonant pump pulse has been devel- 
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FIG. 3: Stimulated Stokes output as a function of the Raman 
shift from the central frequency of the pump pulse. The differ- 
ent panels show the output spectra for various phase offsets. 
The gain oscillates between the S- and A-modes. 



oped both semi-classically |23j and quantum mechani- 
cally [3 Ell- In the quantum case, the Stokes field in 
the non-impulsive regime is generated by spontaneous 
fluctuations of the Raman polarizability. A number of 
approximations can be made to reveal the underlying 
physics. In the standard treatment, one assumes one- 
dimension, plane-wave pulse propagation in the slowly- 
varying envelope approximation (SVEA), and the calcu- 
lation is performed in the transient regime where damp- 
ing can be neglected. The pulse durations (0.1 to 1 ps) 
imply that the SVEA is valid. The assumption of tran- 
sient scattering may be an oversimplification because the 
pulse lengths approach the decoherence times of the C-H 
vibrations!^. With these simplifications, the interaction 
can be reduced to a set of coupled differential equations 
describing the propagation along z of three fields inside 
the medium |27l |: 
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dx 



(1) 
(2) 
(3) 



In these equations. 
Stokes field, and q 
independent variables x (x z and 
reduced space-time coordinates. 

We have extended this treatment to include two vibra- 
tional Raman modes which, in addition to the original 



El is the pump laser field, ss is the 
is the molecular polarizability. The 
T oc (t — zjc) are the 
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FIG. 4: Energy level diagram with two Raman-active vibra- 
tional modes. The three levels are Raman coupled by the 
pump (ujl) and Stokes (ugi and u!s2) fields. The interaction 
generates coherences between pairs of the three levels (gi, 52, 
and ga). An additional interaction is included by allowing a 
two-photon Raman coupling between the two excited vibra- 
tional states driven by either the applied pump field or the 
generated Stokes fields. 



coupling to the same ground state, are Raman coupled to 
each other through two-photon transitions driven by the 
appHed pump field or the generated Stokes fields. The 
relevant energy levels are shown in Fig. ^ where |el) 
and |e2) are the first excited levels of the two vibrational 
modes. Equations2|through|51are the modified equations 
governing the optical fields and molecular polarizabilities 
in the case of two modes, 
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In these equations, the q^'s are the molecular coherences 
for the Raman-active transitions, the Wi's represent the 
population inversion in each of the two modes, and a 
is a constant that describes the relative strength of the 
additional Raman coupling. Specifically, we include the 
quasi-impulsive pump-pump coupling, where the relative 



phase, 4>L, appears implicitly in the terms containing 
£l£l- Finally, although substantial energy conversion 
between pump and Stokes fields is possible, the popula- 
tions of the excited states in our experiments are always 
much less than the ground state, so terms in the calcula- 
tion involving the off-diagonal coupling between the two 
excited states (53) may be suppressed. 

In order to compare our model to the experimental 
results, we numerically integrate the set of coupled dif- 
ferential equations that describe the propagation of the 
fields inside the medium (see Eqs. ^ through . We 
set the amplitudes of the initial coherences (gi and 172) 
to be small random numbers, which in turn determine 
the initial values of the Stokes phases (j)si and (j)s2- The 
equations are integrated over an interaction region com- 
parable to the experimental conditions. At the end of 
the interaction length we calculate the total energy con- 
tained in each of the Raman fields as we vary the phase 
offset (/)L in our model pump pulse. 

We expect the initial phases of the Stokes fields to 
be random because they buildup from spontaneous scat- 
tering. The stimulated Stokes field emerging from the 
methanol cell is clearly multimode, based on its angular 
divergence and far field appearance. When the excita- 
tion volume corresponds to a large value of the Fresnel 
number (F = A/XgL), such multiple spatial modes are 
excited, leading to a large number of independent SRS 
processes in each laser shot [28ll29| . 

The high number of spatial modes implies multiple sets 
of initial phases for the two different Stokes frequencies on 
every laser pulse. Therefore we perform the integration 
over many trials of random initial amplitudes of the co- 
herences qi and (72 (using a white Gaussian noise about 
zero). The control comes from the £l£*l terms in the 
equations for the temporal derivatives of qi . The physi- 
cal picture is that the pump pulse drives the coherences 
gi, which in turn drive the Stokes modes £31 and £52- 
Changing the value of (j)L in the equations changes which 
mode receives a positive contribution and which receives 
a negative. When the phase offset (f)L goes through 7r/2, 
the role of the two modes switches. 

We define the mode asymmetry as the ratio of the dif- 
ference in energy between the two modes divided by the 
sum of the energy in the two modes. Although for a 
given initial pair (51,52) the mode asymmetry can take 
on any value between ±1, on average one of the modes 
is preferentially excited at each value oi <j)L- Figure El 
shows the effect of mode switching for the results from 
both the experiment and simulation. The mode asymme- 
try is plotted along the vertical axis. We find that both 
the experiment and simulation show a variation in mode 
excitation with similar periods. The simulation results 
have been adjusted to the data using two parameters: 
the pump-pump Raman coupling strength (a) and the 
initial phase offset (p^. The depth of modulation in the 
asymmetry increases as the Raman coupling strength in- 
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FIG. 5: Variation in mode excitation as a function of phase 
shift between the two spectral regions for both experiment 
(solid) and simulation (dashed). The vertical axis plots the 
difference in the mode excitation divided by the sum. The 
simulation results have been adjusted to the data using two 
parameters: the initial phase of the field and the coupling 
strength between the two excited states. 

creases (the plot shows the result with a = 7) . Changing 
the phase offset (pL simply shifts the phase of the mod- 
ulation. The phase offset in the experiment is unknown 
due to pulse propagation effects within the liquid. The 
mode control was relatively insensitive to the other pa- 
rameters of the model, which were set to approximate 
experimental values of the Raman threshold. 

In conclusion, we have demonstrated mode selection in 
stimulated Raman scattering in liquid methanol through 
control of the coherence of the pump laser. Control is 
possible even though the pump laser bandwidth is in- 
sufficient to seed the Stokes waves. Using our feedback 
control results as a starting point, we have constructed 
a simple model of the process that could explain the ef- 
fect. The control in our model comes from the Raman 
coupling of the two excited modes by the pump laser. 
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